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ABSTRACT 

We begin by recalling the isothermal, coUisionless, disc-halo. The disc compo- 
nent is the Mestel disc. Subsequently we introduce spiral arms to such an isother- 
mal disc-halo system that are co-moving in the mean with an axi-symmetric 
background. These correspond to a similar disturbance in the halo, which is 
comprised of spiral structures on cones. The arms are necessarily transient due 
to the differential winding in the disc and their gradual destruction is described. 
Although the spiral potentials are weak compared to the axi-symmetric poten- 
tial the arms are not propagating waves on the background, but rather co-move 
with it. They have an effect disproportionate to their relative magnitude on the 
gas distribution in the disc. The gas accumulates on the outside leading edge 
of the 'stellar' arm and an arm-interarm modulation of up to 100% is possible. 
Compatible isothermal, scale-free, distribution functions are found either exactly 
or approximately for all of the coUisionless components of the disc-halo system. 
Repeated episodes of winding arms can produce an exponential disc. 

Subject headings: galaxies:spiral, galaxies:structure,gravitation,spiraI arms 
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1. Introduction 

Our objective in this paper is to construct an isothermal disc-halo and transient, 
isothermal, disc-halo spiral structure from a mixture of collisionless and gaseous matter. 
We restrict ourselves to an infinitely thin disc immersed in a background halo. The initial 
spiral arms are also infinitely thin in one approximation, where they are discrete. Both 
artefacts may be regarded as the result of 'coarse graining' the actual disc and arm. Thick 
discs merging smoothly into halo structures can be stu died in t he sa me fashion, but these 



have already resulted in the well known Evans models lEvans I (1l993l ). The formulation is 
meant to be gravitationally and dynamically self-consistent to a reasonable approximation 
(the CO- moving circular particle velocity is small compared to the disc rotational velocity). 

We begin by summarizing axi-symmetric, self-similar 'isothermal disc-halosQ. The discs 
do not have th e same prob l ems w ith gravitational equilibrium as do rigidly rotating discs 



and arms (e.g. iHenriksen"! f 2 1 ll ) -similar ity class 0). They may require a co r apatible halo 



i n orde r to r emain stable to linear per turbations ( see lEvansfcRead I (1l998al ) , lEvansfcRead 



(!l998bD and iGoodmanfcEvans I fjl999h ). In any case we do find a compatible halo in this 



paper and together, the halo and the disc, define an axi-symmetric 'isothermal disc-halo 
system'. 

There is an infinitely large class (class a = a/5, a positive real number equal to the 
ratio of spatial (1/5) to temporal (1/a) scales) of self-similar rotating thin discs, all of 
which possess differential rotatio n except class zero (r i gid ro tation). This differential motion 
presents the 'winding problem' (iBinney fcTremaine I ( 20081 )) for non axially symmetric 



structures comprised of the same rotating material, which argument implies that such 
structures can not ultimately be stationary. 

This difficulty, which is co mmon to all discs in differentia l rotat ion, has inspired a 
linear theory of spiral structure (iLin. C.C. fc Shu. F.H. I ( Il966l )) (see iBinney fcTremaine 



(120081) for a description of later developments). This theory derives the structure as a wave 
pattern propagating on the disc material. The pattern is assumed to be more nearly in 
rigid rotation with an angular vel o city p .p. However both simulations f SellwoodI f 201ll )) 



and analysis (jBinney fcTremaine I ( 20081 )) suggest that these waves may also be transient. 



This paper takes a rather different approach. The spiral arms are allowed to be 
material arms. They are normally transient and it is this evolution that we study in the 
non-linear limit. There is one case where the arms maybe in rigid rotation and long-lived, 
but after presenting the possibility we do not develop it further in view of the evidence. 

After summarizing the axi-symmetric disc-halo structure that follows from isothermal 
collisionless matter, and discussing what rigidly rotating material arms would have to 
resemble, we model in detail collisionless spiral arms that are co-moving in the mean with 
the background disc rotational velocity. However these do not avoid the winding problem 
and are consequently transient. By focussing on the evolution of the spiral potential as the 
winding proceeds, we describe the gradual destruction by winding of the initially self-similar 



^ These are collisionless systems with similarity class a = 1. 
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arms. The initial arms are maintained so long as the quantity Vt/r {V is the disc rotational 
velocity) is small. This restricts the lifetime at a given radius and the range of radii over 
which the arm persists at a given time. Time is to be measured from the establishment of 
the spiral structure which origin, either by instability or infall, we do not discuss. 

The spiral disc potential must be associated with a non axially-symmetric potential 
component in the halo. The direction in which the causality operates may not always be 
the same, since the transient disc arms may be stimulated by the decaying orbit of an 
inf ailing object. In any case these spiral components are expected to be small compared to 
the axi-symmetric background potential. 

The various components of the disc-halo system are constructed from a scale-free, 
isothermal distribution of collisionless particles plus scale-free, isothermal gas. Considerable 
discussion is given to the boundary condition on the potential at the disc. By using a 
distribution function approach we bypass solving for the detailed orbits of the arm particles. 
The orbits are nevertheless defined by the characteristics of the corresponding distribution 
function. 

The model presented is not a wave theory since the arms are comprised of a 
separate dis tribution of particles that rot ate s in the mean with the disc velocit y. Recent 
simulations. ( Wada.Baba fc Saitoh 1 (j201ll )). (jKawata.Grand fc Cropperl ( 201ll )). and 



observations ( Foyle et al. I (j201ll )) encourage this point of view. The strength of the spiral 



potential is small compared to the total (disc plus halo) axi-symmetric potential, but it can 
nevertheless have a non-linear effect on the distribution of gas in the disc. 

Although we do not solve for the gas dynamics consistently in this paper, it is likely 
that there is substantial streaming of gas and associated magn etic field through th e 
arms. Such streaming can lead to shocks and hydraulic jumps (jMartosfcCox I ( 19981 )) in 



the gaseous matter. The magnetic field is essential to the full understanding of the gas 
dynamics. 

In the next section 2 we derive the isothermal disc-halo solution in axial symmetry. 
This represents the background for the non axially symmetric, isothermal structure. In 
section 3 we discuss the non axially symmetric disc and halo components, including their 
potential and distribution functions. Section 4 constructs an example of the disc-halo 
system with spiral structure. The final section is reserved for discussion and conclusions. 



2. Axially Symmetric Discs and Halos 
2.1. Discs 

It is convenient to describe the familiar axi-symmetric Mestel disc in a differentially 
rotating reference frame. Without axial symmetry such a trick does not work because of 
the winding problem, but in this axi-symmetric example the self-similar analysis in such a 
frame allows us to establish a certain uniqueness for the self-consistent distribution function 
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(DF). Hence the angular velocity of the locally rotating frame is f2 = V jr. 

We do not normally regard the Mestel disc as being an example of self-similarity, but 
in fact it is an isothermal example. The self-similar surface density is uniquely o = S/(5r) 
(S is constant, 6 may be thought of as an inverse arbitrary radius to appear more explicitly 
below), and this yields the corresponding potential above and on the disc, due to the disc, 
as (in cylindrical coordinates) 

27rG'E . 27rGS, , 

'Pd = — ? — arcsmni-) -\ — mor. (Ij 

r 

We obtain this expression in spherical coordinates by letting z r cos 6 and r ^ r sin 6. 
In either expression V^^^ = above the disc and the boundary condition l/{27cG)dz^d = o" 
(or equivalently —l/{27iGr)dg^d = a) is satisfied at the disc. 

The final equation that determines self-consistency is 

cr = y J ^ dvrdv^, (2) 

where an appropriate two dimensional distribution function (DF) must be found. 

A self-similar disc that depends on a constant velocity V (or equivalently on a const ant 
specific energy Ey) falls into the self-similarity class a = 1 (jCarter fc Henriksen I (jl99ll )). 
That is, the temporal scaling 1/a is equal to the spatial scaling 1/5 so that any velocity is 
not in fact scaled. Thus a constant velocity, or equivalently a constant specific energy, is 
compatible with this class of self-similarity. This constant may be used if necessary, to set 
finite limits to the integration over the distribution function. 

By enforci ng rigorous self-similari t y and a steady state in axial symmetry it can be 
shown (see e.g. iHenriksen fc Widrow I f 19951 ) for similar methods and also later sections f 



this paper) a general form of the DF is 

F = K(C)e-(^^+^"^)/*°), (3) 
where E'^ = T + <l>o5i? — the 'energy' in the locally co-moving frame. A constant 
g(y /(2*o)) j^g^g been absorbed in K{C). The one integral that preserves the self-similarity is 
found to be 

C=(r + t;<^)e-(^+^''^\ (4) 
and T = {yl + 'i^|)/2 in the locally rotating frame. 

One readily finds that E'^ + Vv^ = E^, where E^ is the energy in the inertial frame. 
That is Ed = E'^ + ujrv^, which is familiar as the Jacobi integral. 

We note also from the form of Ed that, if v^j, < V, one can drop the term VV(i) in 
equation ([3]). This requires $o > if the bulk of the particles are to satisfy this condition. 
Subsequently $o becomes $oa, which is due to the halo as well as the disc. Using this 
approximation with K constant the DF has the isothermal form 



(5) 
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but with the energy in the locally rotating frame. The integral over velocities continues 
to give the necessary 1/r surface density, and the DF in terms of E'^ is compatible with a 
mean rotation. 

The function K(C) can be arbitrary (we have taken it to be a constant consistent with 
0" = S/5r in the approximate argument above) , since equation (|2]) when the integral exists, 
will always give a oc 1/r. However in general the integral in this equation will yield for S a 
complicated function of V and any amplitude constant appearing in the function K[C). 
For example one might choose K{C) = Ki In (|C| + K2), whereupon with Ki, K2 fixed and 
positive (with \C\ + K2 > 1) cr may be calcul ated in principle give n V an d $0. Such a DF is 
different from that often used in this context iBinney fcTremaine I ( 12008 ). so that even with 
strict self-similarity there is no absolute uniqueness. 

However this apparent generality is spurious if instead a, $0 and V are all fixed, since 
then an arbitrary choice of K{C) ends by defining the associated amplitude constant in 
an ever more complicated way (other constants such as K2 in the example above may be 
chosen for regularity of the DF). It suffices then to make a choice for K{C) that allows 
the chosen values, although this is clearly not a unique choice. A power law in the form 
K{C) = KdV^/C^ serves this purpose and co rresponds to one intuitively composed from 
the energy and angular momentum integrals ( Binney fcTremaine I ( 20081 )). The constant 



is a new fiducial constant with the dimensions of F, and it transpires subsequently that 
q can be any real number smaller than 1. 

We find thus a strictly self-similar DF for the Mestel disc in the locally rotating frame 
according to (|3]) as 



:i + ^ V *o ) 



X exp(-(^^^)), (6) 



o 



where once again e*^^^/^) is absorbed into K^- 

We may calculate a from equation by integrating over velocities. Since the upper 
and lower limits in Vr may be taken as positive and negative infinity respectively, we see 
that g < 1 for the integral to converge. 

The lower limit in bears some thought. A DF of the form (0(x) is the Heaviside 
function) F = Pe~^^ where P = P{vr,Vfj,)Q{V(i, + V) continues to satisfy the self-similar 
Boltzmann equation everywhere, but the same expression without V does not. This means 
that the lower limit should be = —V, when the angular momentum of this particle is 
zero. This implies that C > for all particles in the ensemble. 

The integration over velocities for a now yields (r(x) is the 'gamma' or factorial 
function) 

" - (T^ n'-^)eW^ (7) 
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where 

'l-q)V' 



A 



2$, 



We recall that $o = (27rGS/(5). If all particles were at rest in the rotating frame then 
for equilibrium $o = V^^, if the disc is isolated. However this is not the case for coUisionless 
particles as we calculate below. 

Equation ([7j) is a relation between K^, $o and a plus q. To find the meaning of q 

we calculate some mean quantities. One finds using the DF ([6]) that = $o/(l — q), or 
equivalently 

g = l-=, (9) 

and so the radial dispersion is greater than or less than $o according as g > or g < 0. 
The mean azimuthal velocity is easily found in the same way as 



If the mean velocity is taken to be zero to enforce the net rotation, then setting the 
right-hand side of this last equation equal to V yields the relation between V , $o and q. One 
finds that q < 1 for a reasonable result. To imitate the phenomenon of 'asymmetric drift' 
( Binney fcTremaine I ( 20081 ) ) . one would have to allow q and hence to vary appropriately 



with radius. 

The (squared) azimuthal velocity dispersion is simple in the inertial frame, taking the 
value {v^ + Vy = $o- From this result and equation (fTO!) we find the squared dispersion in 
the locally rotating frame as 

"5-^Ml + ||-2^^|. (11) 



These results parallel those in ( Binney fcTremaine I ( 20081 )). but we have shown that 



the DF may be found by requiring it to be scale-free rather than being an 'ad hoc' function 
of energy and angular momentum. Moreover it follows from this treatment that although 
the DF of the self-similar Mestel disc is not unique, there is a restricted family of possible 
DF's. Each member of this family would give slightly different particle mean dynamics. 
We have chosen a sufficient DF that allows ready calculation and coincides with a previous 
choice. 

The effect of an isothermal spherical dark halo becomes clear after combining the 
argument above with that of the next section. Such a cqmbiri ation was studied long 
ago (e.g. iMonet. Richstone fc Schechter I (jl98ll ) lToomre I ( 19821 )). but we have derived it 



independently from our formal self-similar considerations. 
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2.2. The Disc-Halo of Isothermal Self-Similar Class 

We do not expect the halo above the disc to be in rotation, at least not with the 
amplitude of the disc rotation. We work on the halo therefore in this section in the inertial 
frame assuming any halo rotation to be small. 

A halo that is compatible with the Mestel disc will fall into the same self-similar class, 
but in three spatial dimensions. In three dimensions the velocity and the potential of this 
class are not scaled by powers of r just as is the case for the disc, but the DF is scaled 
slightly differently as 

/ = P(0,^,t;)e-2^^. (12) 
The density scales similarly according to 

p = e(0,^)e-2^^. (13) 

Although the potential does not scale by a power of r, it is always possible when there is a 
constant velocity to include a logarithmic term so that the most general potential may be 
written as 

$ = <l>o(5i? + ^((/),e). (14) 

Here the constant $0 is compatible with the self-similarity because it requires only a 
constant velocity squared for its dimension, while the logarithm {6R) is dimensionless. We 
use spherical polar coordinates and R is once again the logarithmic radius according to 
Sr = e^^. 

With axial symmetry we can ignore the dependence in the potential, and so the 
Poisson equation becomes 

inG ^ 1 d , d-^ , , , 

A spherically symmetric halo of this class gives an inverse square density law (the 
singular isothermal sphere) according to (47rG6s/(5^) = $os and p = Qs/{Sry. With a 
Mestel disc added to this spherical halo the combined potential is = + ^os^R (the 
disc potential is from equation ([T])), that is 

3>c = — — in or -\ — — in Or 

0^ 

-|- — - — (lnsin6' + arcsm/i(cot 6')). (16) 



We write this more simply as 



5>c = — F — inOr-\ — in Or 

0^ 

+ ^![^(ln(l + cos^)). (17) 
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This combined potential of the isothermal disc-halo satisfies the Poisson equation f|T5|) 
if $0 has the combined value 

^oc = {2tiG/5){T. + 2QJ5), (18) 

whence the potential at = 7r/2 is $ocln5r. Normally in disc galaxies the second term 
in the potential is much larger than the first at an appropriate 5, which justifies taking a 
spherically symmetric halo as a first approximation in this expression. 

The question arises as to what type of matter forms the halo? It might be isothermal 
gas or collisionless 'isothermal' matter. However for galaxies it is of some interest to 
consider the constraints that follow from regarding it as comprised, at least in part, of 
collisionless matter. Since the halo density will now depend on the potential through the 
DF, a more complicated dependence on 9 by the potential may be expected. 

We have analyzed the general collisionless Boltzmann equation in spherical symmetry 
in the inertial frame. The rigorous application of a = 1 self-similarity yields the unique 
'isothermal' distribution function as 

P = K^e-2^'^/*°, (19) 

where £h = 'if{4>,9) + {vY/2, and -u is a three vector. When the scaling is applied to obtain 
the physical DF this becomes 

/ = i^,e-2^'^/*°, (20) 

where Eh = \&(<^.^) + ^o^R + (^)^/2. The DF retains this form even in the absence of any 
particular geometric symmetry. 

The density of such collisionless matter is given by 

eh = j Pdvr dve dv^ = (7r$,)3/2A^^e-2*/*°, (21) 

which must form at least part of the density appearing in the Poisson equation for the halo 
(flSj) for self-consistency. In axial symmetry this latter equation becomes ($0 ^ ^oa) 

^^^-(9, + (7r<|.,,)-V2ir,,e-2*'"/*-) 



52 



1 d , ^d^n^ , , 



where 6p represents the coUisional, gaseous isothermal matter in the halo. 

Isothermal gas in static equilibrium satisfies = Ggde~*^2'i'<i/*o)^ where the isothermal 
sound speed must be = $oa/2 to be consistent with self-similarity. The constant 9^^ is 
the value of 0^ at the disc if we take \E'a = there. 



^By 'rigorous' we mean that we have disallowed any R dependence in P which leads to 
the unique self-similar isothermal DF. 
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It is possible to solve equation fl22|) exactly for an axially symmetric disc-halo potential. 
We must use a disc boundary condition 

(23) 



27rG d9 S ' 

In addition we are free to set \E'Q(7r/2) = 0. This allows the isothermal, scale-free, disc-halo 
system to be treated exactly. 

To obtain the solution for ^E'a(^) we introduce the 'ad hoc' constant 

+ ,24) 
and then y = "if a/ ^oa in order to write equation (12^ as 

The solution follows by defining y = u + In sin 6 since the resulting equation readily 
integrates for u{6). We use the boundary conditions u{tt/2) = and we impose the disc 
by {du/dO).„i2 = — 27rGEa/(5$oa). The solution that results for is (a sign ambiguity is 
resolved by requiring to be positive above the disc) 

e2*«/*- = -^sin2ex 
Q + S 

cosh^(c.-VOTSln(^)). (26) 



We have set 



and 



S = (27) 



oa 



cosh^ (C2) = (28) 



As 6^ equation J29D gives (^a/$oa ~ ((1 - y/Q~+^))\ne. This goes to zero so 
that {dg'^)o = only ifBQ + 5 = 1. Otherwise, provided that Q + S > 1, the potential 
goes to positive infinity, the halo density goes to zero according to equation (12T|) . and 
(9e\E'a)o — > —00 as —1/9. This would require a negative mass per unit length on the axis 
and is unphysical by itself. The opposite case when Q + S < 1 has the potential going to 
negative infinity on the axis and the density going to positive infinity there. The mass per 
unit length is then positive, which is also unphysical by itself. 



■^I am obliged to the referee of a previous paper for calling my attention to this. 
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When S + Q = 1 the expression fl26|) simphfies substantially to 



ln(l + v^cos^), (29) 



oa 



where we have added the subscript 'a' throughout this discussion to emphasize that this is 
an axi-symmetric disc-halo potential. This differs from the disc potential included in ( ITTl) 
only by the presence of S, which is indeed equal to unity for an isolated disc. The condition 
S' + Q = 1, is from the various definitions, a useful relation between the gas density plus 
coUisionless density of the halo measured at the disc and the disc surface density. For 
negligible halo gas density this relation becomes 

According to this relation the assumption of isothermal self-similarity in the disc and halo 
implies a kind of disc-halo 'conspiracy', that is a smooth dominance transition in radius. 

Examples of this behaviour of the potential are shown in figure ([T]) when the disc and 
the isothermal gas are moderate perturbations to the isothermal coUisionless halo. We see 
by considering the form for the density that as S decreases, the density contours become 
more spherical. 

This solution has been derived he re in the context of our self-similar is ot hermal 
analysis, but it was known previously ( iMonet. Richstone fc Schechter I f 198ll )). (jToomre 



( 19821 )). although without the isothermal halo gas contribution. 



The distribution function of coUisionless matter at the disc is now comprised of two 
components. From equations (j6]) and ( l20l) we have in fact {Sd{x) is the Dirac function) 



X exp((.-l)^^^^). (31) 



In this expression E^a = {v^ + ^v^ + Vy +vj)/2 + ^oa^R + '^a{(^) for the halo population, 
and E'^^ = {yl + f^)/2 + ^oa^R - /2 for the disc population recalling that = 0. The 
velocities are rel ative to the lo cally rotating frame. The result is similar to the limiting DF 
found by Evans ( Evans I ( 19931 )) when g < 0. 



Such a DF allows for various asymmetries in the stellar velocity ellipsoid at the disc. 
We note in particular that each of Eha and E'^ + Vv^ are constants on their respective 
characteristic. If we follow a joint characteristic (traced by a fictitious particle) by holding 
the DF of equation (|3T1) constant then, provided that v^p/V is small and/or g — )■ 0, we 
should expect each of these quantities to be constant on the joint characteristic defined by / 
constant. Then taking the difference Eha — {E'^ + Vv^) yields that the energy perpendicular 
to the disc 



-j- + ^!a{0) ^ constant, (32) 
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theta 

Fig. 1. — The figure shows the variation of ^ in units of $o in the range of [0, 7r/2] radians. 
The parameters starting from the top curve are S = 0.1, .05, .01 respectively and Q = 1 — S 
in each case. These conditions ensure a relatively unimportant disc in terms of mass and a 
relatively small component of isothermal gas in the halo. 
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along a fictitious characteristi cs close to that of th e thin disc. This appears as a 'third 
integral' in some disc models iBinney fcTremaine I ( l2008l ). In order for the fictions 



characteristic to be close to that of the disc, the disc population should dominate the halo 
population at the disc. The argument is even more direct if we neglect v^/V in the disc DF 
and take g = , so that the isothermal form F = ii'dae"'^^" applies also in the rotating disc. 

This isothermal disc-halo model is theoretically satisfying, but it suffers from two 
conflicts with observations. The first problem is that the strict self- similarity requires the 
disc to be infinite. However, a Mestel disc truncated at a radius Rm has the potential 

^ 27rG'S , r 

*— '"si;: '''' 

to first order in r/Rm (the error is 0{r / RmY). Thus, to this order, the Mestel disc radial 
acceleration shares with that of the spherical halo the property of depending only on the 
mass internal to radius r. The isothermal disc-halo system can thus be regarded as 
forming the central part of some much more extended system. 

The more serious clash with observations is that the mass surface density declines in 
radius as a power law and i iot as an expone ntial. The observed decline of the disc light is 
exponential on large scales f Freeman~l ( 1970[ )) so that one must imagine much dark matter 



in the disc if the Mestel disc is to be taken seriously. This is not generally accepted for the 
following reasons. 

An isothermal HI disc that is supported by a mixture of rotation and pressure has a 
surface density that varies as Tjg{5r)^^ -'^oa)/c^^ obtain the self-similar r dependence we 
must have therefore = + c^. If such a compatible disc were suffi ciently massive, 
i t cou ld provide the dark matter. However a typical HI surface density (iBroeilsfcRhee 



f 19971 )) is IM^/p c^ which, if the sola r neigh bourhood is typical, is a negligible fraction 



of the disc mass fjBinnev fcTremaine I f!2nn8D ). The same conclusion applies to the total 



interstellar medium. 



There is moreover no dynamical evidence for dark matter in our galactic disc ( lEinasto 



f 201ll )). so even a population of low mass stars and remnants (such as black holes and 



neutron stars) is excluded. In the end we are left again with the (steady/axi-symmetric) 
isothermal disc-halo system being physically relevant only within one or two galactic scale 
lengths. Over this range it is possible that the discrepancy between the exponential and 
the power law is difficult to det ect observational ly. This is likely to be particularly true for 
Freeman type II spiral galaxies (iFreeman I (1l970l )). The isothermal halo by itself appears to 
be more widely applicable. 

We do find one promising result in this regard in the succeeding sections, namely that 
wound up transient spiral structure acquires an oscillating exponential behaviour in the 
surface density. This takes the averaged form oc exp \/?)Vt/r for a two-armed spiral at fixed 



'^This is not strictly true for the halo in the presence of the thin disc since then \E'a = '^a{0)i 
but we may regard this as a small effect when a/ ^oa is small. 
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t. Unfortunately this has a rather different shape from a pure exponential in radius, being 
'cuspier', and it only applies to the the spiral structure. Such structure is amplified through 
its effects on the gas however. 

We turn in the next section to study the transient spiral structures that may be 
imposed on an isothermal, axially symmetric, disc-halo background as reviewed above. 



3. Non Axially Symmetric Isothermal Disc-Halo 



3.1. Steady, Rigidly Rotating Structure in the disc 



It may be that recurrent transient spiral structure in galaxies is the rule (e.g. iSellwood 



( 201ll )). and we shall study the non-linear temp oral evolut i on of such arms in the next 



sub-section. However some recent studies (e.g. ISellwoodl (l2012l )) suggest that spiral 
structure may at least occasionally result from growing instability to internal fluctuation. 
In such a self-excited, persistent mode, the resulting arm should have a constant pattern 
angular speed flp if it is to be long-lived when measured in galactic rotation periods. 

Such an arm will be a growing density wave in the background disc, but one expects 
the gravitational influence of an eventual non-linear wave to modify the DF of the stellar 
disc. Thus even the nat ure of the Lindbl ad resonances that are so present in the linear 
theory may be mo dified fISellwood I fl2012D l One way to describe the ultimate non-linear 



development of this process, is to assume the arm to be comprised of particles that have 
been entrained by the wave and move collectively with the constant pattern angular speed. 

The axi-symmetric Kalnajs disc (e.g. iBinney fcTremaine I (j2008l )) is uniformly rotating, 
finite, and has unstable spiral modes especially when rapidly rotating. This suggests the 
constant generation of spiral density waves, but these are likely to be transient rather than 
steady. We would need non-linear evolution into a rigidly rotating material spiral wave, in 
order to have long-lived structure. 

Because of its asymptotic nature and also because of the uniqueness it affords, we 
might assume that the DF of the entrained particles is compatible with rigidly rotating 
self-similarity. The similarity class of a s piral arm r o tating with a non-zero, constant, Qp is 
a = rather than the isothermal a = 1 iHenriksen"! (j201ll ). The scaling of the DF is the 



same in each class, namely F = Pe ^^^^ but the scaling of the surface density a = T,e^^^\ 
The velocities are scaled in the a = class according to -y = Ye^^ with the consequent 
scaling of energy and potential. The logarithmic radius R i s unchanged frorn previous 
sections. By working in the rotating frame it was shown in iHenriksenl f 201ll ) that the 
self-similar DF for a thin disc could be put in the form 

Once again the particle energy in the rotating frame is E'^ = [v^ + f|)/2 + $e//; which is 
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an integral of the particle motion. The effective potential is $e// = $ — Vt^r"^ /2. 

For strict self-similarity appropriate to rigid rotation, one must take K constant and 
the upper limit in energy space Eo-, either zero or oc VL^r"^. The scaled DF P = PQ{£o — S) 
(G is the Heaviside function) remains a solution of the CBE since the scaled energies are 
independent of R. The potential $ is also proportional to in that case, and the surface 
density has the rather singular profile oc r rather than oc r~^. 

The problem with such a material wave is that it must be strictly limited in radius 
and that the halo potential must adjust to stabilize the material with the self-similar form 
( Henriksen I (l201l[ l). The compatible halo is a core of uniform density as perturbed by the 
disc. In the presence of a background isothermal disc-halo the global self-similar 'conspiracy' 
would be broken. This does not seem like a successful model for large scale galactic spiral 
structure although it might describe nuclear structure. We continue to explore transient 
arms in this paper. 



3.2. Transient, Corotating, Spiral Structure in the Disc 

In this section we construct a non-axially symmetric, isothermal structure, that rotates 
with a constant mean circular speed = V. We know that this can not be a steady 
configuration because spiral structure winds up in time due to the differential rotation. For 
this reason we treat the time dependence explicitly. We seek an approximate transient DF 
for the arms plus the details of how it is destroyed in time. 

Such a model conceives the spiral structure to be 'co-moving' with the axi-symmetric, 
inter-arm disc, rather than existing as a linear wave moving on the background. 
This description seems to c orresp ond to the results of recent simulations r eport ed in 
Kawata.Grand fc Cropper"! (1201 ll ) and especially in IWada.Baba fc Saitoh I ( 201ll ). The 



arms (both gaseous and stellar) found in these papers do mainly co-move during their 
transient existence. We do not suggest that all arms behave in this fashion, since we know 
that sufficiently small disturbances will propagate as waves on a background. In fact this 
model might be considered as a non-linear wave, since it does 'propagate' eventually due to 
winding (see figure |2]). Moreover there is likely to be relative motion between the interarm 
gas and these arms while they persist. 

We treat this problem by remaining close to a self-similar evolution in time, at least 
before major winding has occurred. The explicit CBE equation is the disc version of the 
equation studied later for the rotating component of the halo in the next sub-section. We 
have chosen a local frame that is time independent and coincides with the velocity V of 
the flat rotation curve of the background disc. Thus once again Q = V/r and the relevant 
equation becomes 



dtF +VrdrF + (^ - tVrdrVL)d^F 
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y2 

r 

— -(— A / 7 7 —I— 1 1 !' t J .\ r \ I 7^^ 



2nvr + Vrrdrn)d^.F = 0. (35) 



The formal procedure has been discussed elsewhere ( iLe Delliou. Henriksen &: MacMillan 



(120111 ). and references therein) so we will only outline it here. We use a logarithmic time T 
as the self-similar Lie parameter and introduce on dimensional grounds the scaled quantities 
R, Y, ^, \1/ and P according to 

St = e^^, r = Re^^, ^ = (j) + eT,a = J^e'^^ 
F = P{R,^,Vr.v^-T)e-'^,v = Y, 

= $S In (5/2/\/) + $f>T + 

vl/,,(/?,e,;r). (36) 

Formally 5 has the dimension of reciprocal time, but in fact all temporal and spatial 
quantities (and consequently velocities) may be thought of as numerical values in terms of 
some fiducial radius r^, and fiducial time to- 

The form of the potential is equivalent to 

^dr = $S In {6r/V) + ^dr{R, 0), (37) 

and we recall that there is self-similarity in time only if P, E and "^dr are independent of T. 
The winding term destroys this in a secular manner that we discuss below. 

After writing the CBE in terms of these variables we obtain from it in the usual way 
the characteristic equations 

dT ' dT 

dT - '^R^ym)^' ^^^^ 

dYn _ 2VY, V' 

dT R^ R ^ R R ^ 

^ = ~{Y^Yn + VYn + d^^dr). 

We use Y to distinguish the scaled equations, but it is identical to v. 

The characteristic equation may be combined with the characteristic expression for 
dR/dT to give 

A(l„((y,+ me-))^-^^^a,**. (39) 
which in physical variables is the angular momentum equation 

j^{r{v^ + V)) = -d^^ = -d^^dr. (40) 
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The R, ^ characteristics may be combined with the Yr, characteristics to obtain an 
energy equation in the form 



dE' 



dr 



{dr'^dr - SRdR^ildr + e^c^ 



dr) 



+ V{V + v^)(4^{\nRe'^) 



dT 

^dT 

Here E'^^. = ¥^2 + where ^dr = $S In 6R/V + 6^^J^T + ^^r, is the energy in the 
locally rotating frame at the disc. 

We may eliminate Yr/R = d(\n Re^'^) / dT between this energy equation and equation 
(!39|) to obtain 



dT 



dr) 



We note that Edr = E'^^ + ViV + Y^), which is the energy equal to ^dr + [Y^ + Vf/2 + Y^/2 
in the inertial frame but for a constant —V'^/2. Thus the last equation can be written as 
dEdr/dT equal to the terms on the right that involve "^dr- So long as the spiral structure 
remains self-similar, it will be steady in the locally rotating frame. We would like E'^^ to be 
an integral of the particle motion during this phase, and this will be approximately the case 
ifY^<V and if Edr is an integral. 



To obtain Edr as an integral we must set the right-hand side of equation ( H2i) involving 
"^dr to zero. When used with the radial characteristic to eliminate Yr, this is a linear 
equation for "^dr- The general solution has the form 

^dr = ^dr{^-eT + V/{SR),Re'^) 

= ^,,(0 + fi(r)t,r), (43) 

where (f) + fl{r)t = (f)] and (f)j is the inertial frame angle. This merely confirms that a 
steady potential is required to obtain a steady distribution in the inertial frame. An explicit 
dependence on T, or equivalently r in the above expression breaks the self-similarity. 

However we do not wish to describe material in the inertial frame. We can create a 
potential based on transient logarithmic spiral by taking one variable to be a combination 
K, of the above coordinates in the form 

K = ^ - eT + ^ + ( In (fie'^) 

= 0+ (^)lnr + il(r)t. (44) 
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Here we have supposed that radii are in terms of a fiducial quantity that might be Tq = V j 5, 
and we recall that ^ = V jr. 

We retain the other variable as r = Re^^ . Hence 

^dr = ^dr{^^,r), (45) 

and this form must ultimately be made compatible with the Poisson equation. We observe 
once again that although the winding term V/{5R) = Q{r)t is compatible with the 
self-similarity (requires no T dependence), the dependence on r = Re^^ is not. We shall 
see below that this dependence on r is generally required in order to satisfy the Poisson 
equation in the presence of the winding term. Thus the similarity is broken by this effect, 
as might be expected. 

We wish now to write a DF for material at rest on average in the rotating frame, which 
is compatible with the rotating potential. We do not strictly have an integral in this frame 
since equation (H2l) is currently exact in the form 

dE'dr ^ d(V + Y^ 

dT dT ' ^ ' 

which gives the integral E^r = E'^^ + V{V + Y^) as discussed above. However it is clear that 
provided Y^/V in the local rotating frame is small, we may treat E'^^ as an integral for those 
particles. This condition generally holds for the majority of particles in spiral galaxies. 

Since this is the only identified integral, we write the characteristic solution of the 
Boltzmann equation as the approximate DF P = F{E'^^)e^'^ . But H = J P dYjidY^ and 
this should be independent of T for self-similarity, to which behaviour we wish to remain 
as close as possible for uniqueness. Thus, recalling the form of the potential (e.g. see after 
equation fj4T]) ). we see that we must have the isothermal DF in the locally rotating frame 

F{E'J = Kdre^p(^-^y (47) 

where K^r is the normalization for the transient spiral distribution function. In order for 
the bulk of the particles to obey the condition Y^ < V we should require $fj < V\ The 
mean velocity of these particles is zero in the locally rotating frame, due to the symmetry 
of the DF. When other components are present that may be described by the isothermal 
DF, the potential in the exponential will be the sum of the various potentials. 

This DF also gives S oc 1/i? and hence /r as it should for self-similarity. 

However because of the likely dependence on r in "^dr that we now pursue, this self-similar 
behaviour will be broken in general as the winding continues. This also leads to a more 
interesting radial variation of the spiral surface density, although it is transient. 

To obtain the equation for the disc potential above the plane we use the Poisson 
equation 

^ -de{smede^dr) + -^9l^dr) = 0, (48) 



sin 6 sin 9 
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and insert the spiral form 
to find eventually 



= $i;hnr + *d,(«;,^,r), (49) 



\ r sm y J 

+ ^a,(sine9,*d,) = 0. (50) 
sm d 

For brevity subsequently we write the differential operator in this equation according to 

+ $S = 0, (51) 

and may be taken zero by absorbing the log potential into the corresponding 
axi-symmetric term, when present. 

One can only neglect the r dependence in this equation, and so preserve strict 
self-similarity, if Vtjr = fl{r)t < e/6. This might have been expected, but the dependence 
on the initial winding angle e/S is of interest. In the circular arm initial limit {e/S — > oo) 
the winding can, not surprisingly, be for an indefinite time. 

This linear equation is readily solved in terms of modes of the form 

= In sin ^ + e^'"'^^T{e)n{r), (52) 

where as already remarked may be absorbed into and so taken zero here. We 
recall that k = (p + {e/d)\nr -\- Vt/r and is in the locally rotating frame. The log spiral 
is completely wound up at a fixed r when fl{r)t = 27r. This gives ~ 10^'^ years at r = 10 
kpc and V — 200 km/sec. However there is an outward moving 'winding wave', given by 
Vt/r = est < e/S < 27r, outside of which the log spiral remains recognizable and similarity 
is maintained. 

The rotating DF and the rotating potential are linked through the disc boundary 
condition ^ 

27rGadr = — {d0'^dr)\e=n/2, (53) 



where 



adr = I F{E',) dYndY^ = 27r$2^*^— • (54) 



e do 



This condition is generally difficult to satisfy for all k, for a single mode, because of 
the exponential dependence of a on the potential that follows from the last integral. 
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Fortunately we can choose to satisfy it everywhere by adding isothermal gas to the disc that 
is not described by the isothermal collisionless DF. If however the arms are deemed to be 
predominantly comprised of collisionles matter, then we can satisfy the boundary condition 
only at discrete values of k. These then become the idealized spiral arms much in the same 
fashion that the razor-thin disc is an idealized thick disc. 

For a fixed spiral trajectory, the particles ought to be constrained to move 
one dimensionally along the spiral in the co-moving frame. However at fixed /t = 
+ (e/5) ln(r) + Vt/r = v + Vt/r, the log spiral sweeps over different trajectories according 
to dv = —Vd{t/r). The effective arm therefore will be extended and distorted so we 
continue to allow a two-dimensional DF to describe the arm particles. This is illustrated 
in figure (|2]). We see that the winding occurs early at small radii and later at larger radii. 
Moreover the arm appears to move as a non-linear wave in the co-moving frame at large 
radii before it is completely distorted, which one expects to thicken the actual arm. The 
slowest destruction occurs for the larger winding angle so that arms with large winding 
angles (including rings) are most likely to be observed. 

This concludes our model for transient arms in an isothermal, thin, disc. In the 



Fig. 2. — The figure on the left shows the k = spiral with initial winding angle e/5 = 1 at 
different times in the xy plane. The figure on the right shows the same spiral in the xy plane 
with initial winding angle e/6 = 0.5 at the same times. The times are from the top curve to 
the bottom at x = 1; Vt = 0,0.05,0.1,0.3,0.5,0.7,1.0 respectively. The initial spiral with 
smaller winding angle is distorted more rapidly. Distances are measured in terms of some 
fiducial radius Tq. 




0^ 
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next sub-sections we treat the compatible non axially symmetric halo and the necessary 
isothermal gas component. 



3.3. The Transient Non-Axially Symmetric Halo 

The total potential in a disc-halo system must satisfy the Poisson equation in the form 

V^<|. = AnGp, + 4vrGS^^^i^^^, (55) 

r 

where the total potential is the sum of that due to the disc and that due to the halo namely 
$ = + $/^. This equation separates for the two components to give 

V2$, = 0, (56) 

where the total disc potential satisfies the boundary condition (155]) when a is the total 
surface density. 

Recalling the previous sections, we can form the potential due to the disc from two 
components. These are respectively the axi-symmetric potential $da(^) given as the disc 

component in f|T7|) plus the rotating structure = Inr + "^dr- The potential "^dr is a 



solution of equation (I3U]) . Thus we write {5 of the steady state spatial scaling is replaced by 
the time dependent scaling S/V) 

$d = — r-— (lnr + ln(l + cos6')) 

0/1/ 

+ $i;hnr + ^,,(K,^,r), (57) 

where is the axi-symmetric surface density so that the first two terms in this potential 
comprise ^da- We will use for the rotating, non-axi-symmetric component. Both E^^ 
and Y^dr will include an isothermal gaseous component in general, and so indeed may 
the volume density p^. We adopt the notation ^do = ^do + {^irGEdaV/d) = + 
subsequently. We note that each of ^da and ^dr satisfy a Laplace equation above the plane. 

To the extent that the spiral structure is a small component of the disc-halo, we 
might expect the first approximation for the halo potential to be an axisymmetric function 

^ha{0)- This may be written as ^ha = ^^ho^"^^ + ^/ia(6'). The combined solution for 
= "^hai^) + ^da(6') is givcu by equation (1291) . Hence "^fhai^) may be found by subtraction 

as 

^ha = - ^da (58) 

= $oJn(l + v^cos^) -$i^Mn(l + cos^) 
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We recall that S = ($^"''/'^oa)^ in current notation, so that dg'^ha\-K/2 = as it should. The 



factor = (^^fl + $1"'' in current notation. 



ha\-K/2 

do ' ho 

However the rotating spiral disc structure will impose a rotating perturbation on the 
halo potential/density (^hr{i^iO,r) = 6*, r) + ^j^'^lnr, which form is derived below. 

It may be that the the causal order is inverted. This would mean that the 'spiral' halo 
structure is actually the origin of the disc spiral structure, being itself due for example to 
the decaying orbit of a merging object. We do not have to decide this point here, as the 
formalism is the same in either event. 

Consequently we write for the halo potential $/i = ^ha + ^hr 

$ft = ^ho Inr + ^hr{K, e, r) + ^^taiO), (59) 

where ^ho = ^h'o + ^ho- "^^^ Poisson equation for the total potential becomes (after adding 
the Laplace equations for and and considering 9 < n/2) 

v^i^Ue) + ^da{e)) + v^{<^dr{K, e, r) + $;,,(/€, e, r)) 

Here Q is as in equation except that 6 ^ S/V, and $o — t- $o„ = ^^o + ^ho- 

The axisymmetric part of this last equation (using $„ = + ^ha and in Q) 
satisfies 

V ($a) = — ^ e (61) 

and has the solution fl29p for when the disc boundary condition is imposed. As remarked 
above this solution is the simplest approximation to the halo potential, which follows by 
setting the rotating potential components of the disc and halo equal to zero. The next 
approximation may be found by expanding the second exponential in equation (160|1 . by 

neglecting = ^to + ^d^o ^oa, and by subsequently using equation fl6T]) . to find the 
inhomogeneous linear equation 

V\^hr + ^dr) = -2^e<^){^hr + ^dr). (62) 

In this equation "ifdr is known from equation (150|) . together with the disc boundary condition 
in terms of T^dr- The exponential is known from equation (129|) with Q + S = 1. 

Our task is now to describe the coUisionless material comprising the rotating isothermal 
halo component. This argument parallels our discussion for the rotating disc component 
in the previous section, but must be done in spherical geometry with time dependence. 
Although it is cumbersome, we state here the complete CBE for such a problem. The 
reduction to the disc CBE is immediate by setting the 3D DF / = F6d{0 — Tr/2)6D{vg) {6d 
is the Dirac function) and integating from 7r/2 — eto7r/2 + e over theta and from — e to e 
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over fe, and letting e — )■ 0. The equation is (recall that Vl = V/r) in the locally rotating 
frame 

dtf + Vrdrf + —def + ( — ^ - Vrtdrd] d^f 
r Vrsmfc' / 

+ 2fi sin (e)t;^ + n\ sin^ (0) - 9,$;,,^ 9,,,/ 
cot ^ + 2nv^ cos ^ + sin (^) cos (^) - ^ - -de^hr 

2nve cos ^ + ^ cot ^ + ^ + nvr sin ^ + -^—d^^hr 
r r r sin 6* 

X d,J = 0. (63) 

We convert this equation to self-similar variables in the usual way by assigning (the T 
dependence is because of the gradual destruction of the self-similarity by winding) 

6t = e^^, r = Re^^, ^ = (j) + eT,v = Y, 
f = P{R,^,9,Y-T)e-''^, 

= 

^tI{'^nR + 6T) + ^UR,^,0;T), 

p = Q{R,^,9-T)e-^^^,Q= [ P d^Y. (64) 



The equation that results from inserting these variables into equation ( l63l) has the 
following characteristics: 

dT ' dT ^ ' dT R' 

dT RsinO 6RR' 

dYn Yn^ + Y? 2V 

-T^ = — ^ + — F.sin^ + — sin2^-9Ar, 

dT 2 R ^ R " ' 

dY Y^ 2V V"^ 

^ = ^cot9 + —Y^cose + —sm{9)cos{9) 

al K K K 

- ^-i9e<^^r, (65) 
dY,, 2VY, „ Y,y, „ YrYs 



cos 9 cot 9 — 



dT R R R 
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We can combine the R, 9 and characteristics to obtain 

^(/2sin {9)Y^ + VRsm\e)) = 

- dR{Y^sme + VsmH)-d^^hr, (66) 
which, by returning to physical coordinates becomes the angular momentum equation 

-^{rsmie)v^ + Vrsin' (6)) = -d^^hr- (67) 

Just as is the case for the disc we can combine these characteristics to obtain a relation 
for the change in the inertial frame energy along a trajectory as 

di 

( V d\nR\ ^ ^ 

Here 

Ehr = El + V sin (6) {Y^ + ^ sin 6) , (69) 
where the energy in the locally co-moving frame is 



El = ^ + ^hr. (70) 



We substitute the form of ^hr from equation flMj) into the right-hand side of equation 
f l68|) . and in order to obtain the conservation of inertial energy we set the resulting expression 
to zero. This yields the compatible form of "^hr, namely \l//ir(r, 6,^ — eT + V/ (SR)). As for 
the rotating disc component, we can incorporate the dependence on a transient logarithmic 
spiral by introducing the variable k to write 

<i/hr = '^hr{r,e,n). (71) 

This justifies the form of the rotating halo potential that we used at the beginning of this 
section. An explicit dependence on r = Re^^ destroys the self-similarity. 

The DF for this halo component follows from P = fe^^"^ where / can only depend 
on integral constants. To describe a structure in net rotation we wish it to be a function 
of E^. However this energy is only constant according to equation ( 169|) for Y^ < V and 
d9/dT small. The 'small' must be with respect to so that using the theta characteristic 
we require {d6/dT)/Q = Yg/V<l. Thus for consistency the DF must decline rapidly when 
Yfp and Yg exceed V. 



The isothermal self-similar form is clearly necessary as 

-(%) 

/ = Khre <o , (72) 
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since it succeeds in producing the expected p oc with no exphcit dependence of P 
on T. This self-similarity is broken through the dependence on r = Re^^ in '^hr^ as may 
be required by the Poisson equation. When calculating the collisionless halo density from 

the isothermal DF, the total potential must be used and ^ (^ho + ^do = ^o- By our 
approximations we should have > $„. 

We turn next to the expressions for the various potential components. 



3.4. Non Axially Symmetric Potential solutions 

The basic potential for the disc-halo system is given by equation ( 12^ for the 
axi-symmetric component, and the solution of equation ( l50l) for the spiral disc component. 
After solving for the disc spiral component, the non-axially-symmetric halo component may 
be found in principle from equation ( l60l) . or approximately from equation (1621) . 



However equation ( 1621) is only readily separable using an ansatz of the form ( l52l) . if 

one takes $1'^'' = $1^'^ + ^^J^ = and absorbs the log term into the axi-symmetric potential. 
Equation ( l62l) then takes the form 

~ + V^C0S^)2^" ^^^^ 

where \E'r = "^dr + ^/ir- This equation can be solved in separated form, but it is essentially 
only known as a series. We reserve the complete exploration of such halo spiral structure to 
another work. 

We can only expect to find regular spiral structure in the disc when Vt/r < e/6, so that 
our discussion of regular spiral arms will be restricted to that limit. However to identify a 
possible description of the evolution of the arms we proceed briefly with the general case. 

Equation ( 150|) with the modal ansatz of equation ( 152|) is resolved into two equations 
(e ^ e/6) 



= 0, (74) 

= 0, (75) 

where fc^ is the separation constant (positive or negative in general) and ( = Vt/r. Near 
( = the appropriate solution of the second equation is 7^ = 1 and fc^ = 0. The solution 
to the first equation is then simply found in terms of associated legendre functions and so 

$-(«:, ^)-$fj In sin ^^v^X 



- 26 - 



=e(^'"«)(Ci™ PZXx) + Cm Q™ .(x)), (76) 

where and denote the associated Legendre functions and Cim, are complex 
modal constants. This potential, together with the axially symmetric disc-halo potential, 
will be our principal concern below, but it is of some interest to examine the evolving radial 
dependence. 

The modal solution for the radial dependence takes the form 

7^„(C) = exp^(-^ + ^lnC)(^„^MA,^(y3mC) 

+ A^^WxA^mQ), (77) 



where M, W are Whittaker functions, A = \/^me/2, and /i = A/(r^M?7Ti)M-^A!|^/2. 

We can simplify this expression somewhat by considering the radial dependence near 
C = e, where it should represent a rapidly winding spiral. In this limit the radial equation 
becomes approximately dPTZ/dC,"^ = {k'^/e^ — rn?)lZ. Hence k"^ < m^e^ implies an oscillation 
in C (i.e. 1/r at fixed time) corresponding to a winding of the arm. The value /c^ = e^m? 
corresponds to a marginally stable case, where the deviation from an arm intially constant 
on K = constant is linear in 1/r. In this limit /i = a/ (3me — i){me + i) /2 in the Whittaker 
functions. 

The asymptotic behaviour of the Whittaker functions at large argument are 

_ r(l + 2/x) e^'"^/^ 
^''^ ^ r(l/2 + /i-A) (v/3mC)^' 

Consequently, it is the Whittaker M function that describes the destruction of the spiral 
arm with increasing C,. It is of interest that this destructive evolution produces an oscillating 
exponential decreasing with increasing radius along the arm. The amplitude is proportional 

to (r/m\/t)^e("^*/'^)^/2 

The modal analysis for the function ^1^^ based on equation ( TfSll yields the two equations 
f l7^ [75]) : but with the additional term 2(1 — S)/{1 + v^cos6')^ in the bracket multiplying 
T(6'), in the first of these equations. The resulting equation is solvable formally in terms of 
a Heun series, but it is best studied numerically. A non-trivial exception is when the disc 
dominates the halo so that S* ~ 1. Then the halo spiral structure satisfies the homogeneous 
equation fl50l) . but with different boundary conditions. A sum over modes might be required 
to describe the orbit of an infalling object, but the separated modal form of equations (17^ 
~5]) is relevant if the disc spiral is the origin of the halo disturbance. 



We are now equipped to consider in the next section the properties of transient spiral 
arms embedded in an axi-symmetric, isothermal, disc-halo. 
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4. Transient Spiral Arms in an Isothermal Disc-Halo System 

We study in this section examples of 'initial' (the creation of the disturbance does not 
concern us here) spiral arms in a disc-halo system. Various components are considered. 
These are comprised of coUisionless particles with the corresponding isothermal DF, and/or 
isothermal gas. There is the axi-symmetric disc-halo background, the spiral disturbance in 
the disc associated with the spiral arms themselves and the consequent spiral disturbance in 
the halo. The latter component does not appear in the following disc boundary condition, 
but it could be observable in edge-on galaxies. Its simplest form would be a series of 
logarithmic spirals on cones with amplitude decreasing with decreasing 9. 

CoUisional material must be present, in order to satisfy everywhere the boundary 
condition ^ 

Here aa is the axi-symmetric background density, adr is the rotating spiral density, and 
CTg is an isothermal gas density. We can calculate 0"^ and adr from their corresponding 
distribution functions in terms of the potential components that were presented in the last 
section. The gas density is taken normally to be determined by this boundary condition. 

In fact the gas distribution is subject to the same potentials as are the other disc 
components and in principle its velocity and density are determined by the hydrodynamic 
equations. However the gas behaviour is subject to the magnetic field, particularly in 
the inter-arm regions. Thus, in the absence of major streaming, the gas is likely to be in 
(isothermal) magneto-hydrodynamic (MHD) quasi-equilbrium. By fixing the gas density 
from the boundary condition ( 1791) . we are effectively determining the (quasi, because of the 
winding spiral arms) quasi- equilibrium magnetic field. Such a field in the disc would have 

the equilibrium form B = be~^^, where b = b{n, r). This raises the possibility of comparing 
the consequent magnetic field structure with observations. But we leave this aspect to a 
future work as the required MHD equations are formidable. An eventual complete disc 
solution will require this gap to be closed. 



The boundary condition fl79|) reduces to 

1 

2^Gr 

when we recall that is defined so that 

1 



CTdr + (ygr = - WZTTzi^e"^ dr\n/2) , (80) 



an 



/2J 



2nGr 

Because we take \&a(7r/2) = 0, and because we may choose \E'^r(7r/2) = as well as 
and 



absorbing ^^J^ and ^^f^^ into $00, we can write 



cJd, = 27r$„,— e"^. (82) 
r 

For the gas density we take agr = T,gr{K,r)/{Sr) where for the initial spiral structure the 
dependence on r is ignorable, just as for the potential. 
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We begin our investigation of the boundary condition f lHOj) by considering the solution 
fl76|) for the initial spiral disc potential in more detail. 

We only consider one mode at a time in this treatment, normally m = 2 since this 
correponds to many observed spirals. One might retain a full fourier analysis of in order 
to satisfy the boundary condition flHOl) over a range of k, but one mode would still have to 
be dominant in order to match the observations. 

The constants in the solution f l76|) may be chosen freely, but it is useful to check that 
our symmetry requirement ^E'^rl^) — ^dr(^ ~ ^) satisfied . If we consider \E'2J- to increase 
away from the disc towards the axis, then this symmetry enforces our assumed asymmetric 
boundary condition {do'^]^^)T,/2- = —{dg'^]^j,)T,/2+). Moreover the gravitational acceleration 
of the disc is then towards the disc. 

The Legendre fun c tions at positive x, using the 'cut' employed for example in 



Gradshteyn fc Ryzhik I (Il994l ). have a non trivial relation to those at negative x. To enforce 
the symmetry, and to retain two free constants, the constants Cim and niust be related 
to the constants at x = 0—, namely C{^, Cg^, by the relations 

_ 71 
Clm = C'lmCOS0c - - Sin0c C2m, 

2 

C2m = sin (pc Cim - COS (f)cC2m- (83) 

IT 

We have defined the complex angle 0c = (1 + it)rmT. 

It is possible to insist that the constants retain their values across the disc because 
equations (!83|) then become homogeneous with a zero determinant. However in such a case 
the ratio of the constants is defined in the form 

C2m 2 sin 0c 2i sinh {men) ^^^^ 

Cim TT 1 + COS0C TT (— 1)™ + cosh (mevr) 

Once the winding angle of the spiral disturbance is fixed, this condition reduces the free 
constants to one. Moreover we can suppose that Cim is real since any phase constant will 
simply add an arbitrary phase to e*"^". However this reduced case may not allow us to have 
spiral arms of arbitrary amplitude. 

We have taken the k dependence of our mode to be periodic, as is customary in linear 
wave descriptions. However in a non-linear treatment aperiodic solutions may also be 
possible. These would require m = —ip where p is a real number, and so the potential 
would be aperiodic. This implies discontinuities in the spiral disc potential. In a non-linear 
disturbance after coarse graining such discontinuities may be realized as collisionless 
'shocks'. In resolved detail they would be regions of rapidly changing potential and surface 
density, probably involving normal gas shocks. They might be expected along the edges of 
the spiral arms, but we shall not consider this possibility further in this paper. 

The boundary condition (!80|) becomes explicitly {x = cos 6) 
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= (85) 

where ^I/2r(0) given by the real part of equation ( 1761) and the real part of the 
derivative is also known from this expression. We have reintroduced units of length here so 
that G may have its normal dimensions and value. 

This last equation can be satisfied everywhere by taking it to be an equation for 
S™(fi;), provided that the constants may be chosen so that the net gas density + (S^ 
is re-labelled for consistency) is always positive. 

We take as an example the asymmetric mode m = 2 with an initial winding angle 
e/5 = 2, which corresponds to the complementary pitch angle of 26°. 56. Then the boundary 
condition becomes (dropping the m = 2 superscript) 



^9r(,K) = 7^37777 l^oaVlDP COS (0DP + 0i2 + 2k) 



5_ 

2nGV' 

+ ^oay2DQ COS (0DQ + 022 + 2k)] 

= -2'K^oaKdr exp [-?/lP COS (0p + 0i2 + 2k)] 

X exp[-y2Q cos (0q + 022 + 2k)] , (86) 
where we have introduced the 'ad hoc' definitions 

Plm ^co.W + 2^)f-^^J^^^Pe^^^^ 
yj-n V[2i) 

Qlm ^2v/^sin(7r(l + 2^))^^^^^t^ = ge^*S 
a.(P^)|o ism(.(l + 2.))I|l±|.Z.Pe^^^^ 
d.{QlM V^cos(7r(l + 2^))^||ii^^Dge^*-«, 

in terms of the phases and absolute values of the associated Legendre functions. In addition 
Hi = |Ci2|/$oa and ?/2 = \C22\l^oa- With the moduli and phases of C12, C22 given, plus the 
real values for $oa and Kdr-, equation ( !86|) determines the required non-axially symmetric 
gas density. A simpler example is given by equation fl84|) . which for the present example 
gives C22/C11 = — 2^/77 very nearly, so that with 0i2 = 0, 022 = — 7r/2 and y2/yi = 2/tt. 
This leaves only $oa, Kdr and yi to be assigned. 

This procedure does have the merit of indicating that, without gas, spiral structure is 
not possible in this theory. However there are hidden consequences to be explored, since in 
the axi-symmetric disc-halo the values of Sg^, $0^, and corresponding halo quantities are all 
constrained by 5* + Q = 1- Moreover, this approach makes the dynamics of the gas subject 
to this boundary condition, which dynamics we do not explore further in this work. 

It is however instructive to consider the other extreme, wherein the non-axially 
symmetric gas behaviour is unimportant. This can be done by taking its value in equation 



- 30 - 



f l86|) to be constant. We find then, as was indicated earlier, that the one mode boundary 
condition is only readily satisfied at discrete values of k. These discrete values then 
delineate the 'arms'. Equation fl86|) must then be satisfied with constant in the arms, 
while between the arms we consider only the ipay^ga pair. 

To be consistent we must neglect the effect of dx'^drlo on the surface density in the 
inter-arm region. For example let us suppose that the 'gas' is in pure rotation with a 
constant linear speed V^, and that it obeys the isothermal self-similarity generally. This 
requires ag = T.g{K)e~^^, and the 'pressure' in the surface due to this component to have the 
form ps = T[K)e''^^. Then by applying the radial and azimuthal equations of equilibrium 
to the gas in the inertial frame, we obtain respectively (prime denotes d/dn and we include 
neither viscosity nor magnetic field) 

r'--J = \^giy^-^oa). (87) 

and 

By neglecting '^^r and in addition by requiring = S^^ and hence r to be constant, we 
see that a simple solution gives 



^ + V^- $oa = (89) 

In general the gas dynamics is complex, and if we accept the value of the surface 
density from the boundary condition (l86ll . we would have to solve the gas equations for the 
magnetic field if the gas is in equilibrium. There may well be a fiow field in addition to the 
magnetic field and indeed time dependence, all to be rendered consistent with the boundary 
induced value. But this problem is decoupled from the problem we solve here. 

For definiteness we study the numerical solution of the boundary condition under the 
simplifying assumptions of equation (|8^ . It happens that, for a two-armed mode with a 
large winding angle that we study here, the condition of equation (jHlj) is very nearly the 
same as requiring \l/dr(/^, 0)/$oa = 0. There is a difference in the two conditions at the level 
of a part in 10"^, which allows \l/rfr-(fi;,0) to be small but not strictly zero. We use this 
example here for illustration, but in general '^dr niay not be as small, since the constants 
are arbitrary provided equations ( |83ll are satisfied. 



An inspection of the boundary condition under these conditions confirms that the 
constants may be chosen so as to satisfy it only at discrete values of k. We choose these to 
be at /€ = and at k = vr for the two-armed case. The boundary condition ( |86ll becomes 
thus at K = 



— — + 7 — Kdr exp 



2 

—P COS (j)p Q sin I 

TT 



yi 



2 

yi{DPcos(j)DP + -DQsm(j)DQ), (90) 
vr 
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whichis an equation for yi. This latter measures the strength of the spiral potential relative 
to the axi-symmetric potential. The parameter E = {2TTGV/6)(T,gr/^oa) measures the 
surface density of gas in the spiral arm while the parameter A = {4:7r^GV/S)Kdr is a 
measure of the collisionless surface density in the arm, both essentially as a fraction of the 
axi-symmetric surface density. The bracketed expressions that depend on the amplitudes 
and phases of the associated Legendre functions are pure numbers. 

If we write this last expression schematically as 

E + Ae~^y' = Byi, (91) 

then B = DP cos (pop + ^DQ sin (pDQ ^ 2.48 x 10^ and C = —P cos (f)p — sin (1)q ^ 
—1.6 X 10^. The parts depending on Qli{0) are negligible. 

This schematic relation is easily solved given values for E and A. Relatively large 
choices such as A = 0.5 and E = 0.05 or vice versa give yi ~ —2.2 x 10~^. At A = OA and 
E = 0.01 or vice versa give yi ~ —4.4 x 10"''. Either choice renders the spiral potential 
small in terms of the axi-symmetric structure of the disc as is required for inter-arm 
consistency. 

The opposite extreme, where the gas distribution is essential, can be illustrated by 
solving the boundary condition (186|) for Tjgr{K) under the approximation of equation fjH^ . 
In this procedure we assign yi as well as the parameters E and A. 

Thus in terms of a schematic notation we have, 

E{k) = yi{B cos 2k + Bl sin 2k) 
-Aexp{yi{-Ccos2K + Clsin2K)), (92) 

where Bl = —DP sin (p^p + ■^DQ cos (poQ ~ 7.5 x 10^ and CI = Psin0p — ^Qcos(f)Q ^ 
5.3 X 10^. This yields the function E{yi, A) for the scaled spiral gas density that is required 
ClS db function of the physical parameters. 

In figure ([3]) we show two examples of the gas density dependence on k for A = 0.01 
(upper large amplitude curve), A = 0.1 (lower large amplitude curve) and yi = A x 10~^. 
For larger values of yi the amplitude of the spiral oscillation tends to dominate the 
axi-symmetric density. Larger values of A render the oscillation more asymmetric about 
zero, with the minima coming to dominate the axi-symmetric gas density. 

The smaller pair of curves on the figure show the collisionless particle density variation 
for A = 0.1 (top curve) and A = 0.01 (bottom curve-that has been multiplied by five for 
visibility)) for the same value of yi as for the gas density. We see that the relatively small 
variation in the particle density (and potential) leads to a magnified reaction in the gas 
density variation by more than a factor ten. The gas density peak is slightly leading (larger 
(f)) the particle density peak at a given radius and is slightly outside (larger Inr. This 
reverses as the winding proceeds, as can be seen from figure ([2]) 

This concludes our model for transient spiral arms in the thin disc limit. Many 
variations of the model are possible if C2 is decoupled from Ci. The principal characteristic 
of this model is that the spiral arm is co-moving with the background disc until it is 
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Fig. 3. — The larger pair of curves is the non axially-symmetric gas density as a function 
of K for t/i = 4 X 10~^ with A — 0.01 at top and ^4 = 0.1 underneath. The smaller pair 
of curves give the non axially-symmetric coUisionless particle density for the same value of 
yi and A = 0.1 above and five times the curve for A = 0.01 below. Larger values of yi 
will produce negative total gas density. The gas density responds very strongly to the weak 
spiral potential. 
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destroyed by winding. It is perhaps worth remarking also that an oscillating exponential 
decline appears in the spiral potential as the arm is wound up, which would lead to similar 
transient behaviour in the surface density of the disc. 



4.1. Summary of the Distribution Function 

The various components of our disc-halo system have been described in terms of their 
individual distribution functions. The question arises as to whether the sum of these 
distribution functions is 'valid' (i.e. satisfies the CBE) description of the whole system. 
Fortunately, under the approximation that in the locally rotating frame v^f, < V, the 
axi-symmetric and non axi-symmetric distribution functions take the same isothermal form. 
This holds for the disc and for the halo. Hence writing the appropriate isothermal DF with 
the total potential, remains a solution of the CBE for every component. 

The sum DF is comprised of two halo components and two disc components. If we use 
the approximation wherein the comoving v^p < V (the exact DF is the disc part of equation 
(!3T|) ) these take the form 

E,, e', 

+ KdaC-"^ + KdrC-^, (93) 
where (all velocities are in the comoving frame) 

E'.r = ^ + t(>r/2). 

In these expressions the total potential is 

$ = $oJnr + ^a{0) + ^dr(/t,r) + ^/,.(K,r), (95) 

with the r dependence only developing as the winding continues and \E'a(7r/2) = = 
\E'hr(7r/2). To the extent that "^dr and "^hr are small compared to $oa they might be 
neglected for the coUisionless particles. However we have seen that these small potentials 
can have a major influence on the gas distribution. Moreover they are essential to the spiral 
nature of the disc and halo. 



The distribution function approach that we have used for each component avoids 
the question of the actual particle orbits. These may be found in principle from the 
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characteristic equations of the non axi-symmetric Boltzmann equations for the disc f l38|) 
and for the halo fl65l) . The initial state orbits, before significant winding, can be studied by 
neglecting the terms in Vt/r. The halo spiral density disturbance can be found in principle 
by using (^hr in the Poisson equation 

V'^^hr = ^rrGphr. (96) 

The detailed study of the corresponding orbits must await another work, but it is clear that 
the resulting spiral distortions are of interest as possible infalling orbits. 



5. Discussion and Conclusions 

We have studied the construction of spiral arms and discs based on the distribution 
functions that are dictated largely by isothermal self-similarity. Section (2.1) incorporates 
the axially symmetric Mestel disc into this scaling class, by using the frame with constant 
rotational velocity . This discussion leads to section (2.2) where the compatible (collisionless) 
halo is studied in some detail in the inertial frame. The solution for the isothermal disc-halo 
potential is given in equation ( l29l) and the combined disc-halo distribution function is given 
in equation ( 13T|) . It implies an approximate integral in terms of the energy normal to the 
disc, this solution forms the background for the spiral structure. 

Most of our new results are to be found in section (3). Here we treat a spiral arm that 
is comoving with the isothermal background disc. Consequently it is subject to secular 
winding and hence is transient. By treating isothermal self-similarity in a time dependent 
fashion, we were able to show the effect of the winding on the distribution function and on 
its potential. This winding may be neglected up to a certain time at a certain radius, which 
time increases directly with radius. Thus the transient arm is perturbed from the inside 
out. 

Beyond a critical radius at a given time, the distribution function remains isothermal. 
The potential is required to be a function only of the spiral coordinate initial 
condition, but it becomes progressively dependent on radius as the winding destroys the 
self-similarity of the arm. As the destruction proceeds the non-axially symmetric potential 
adopts an oscillating exponential behaviour. This can be significant in the gas distribution 
by non-linear amplification. If so then after many episodes of transient spirals the isothermal 
disc will become exponential, although not in the Sersic form. 

Another unusual element of the model is the necessary, isothermal non axially 
symmetric structure in the halo. This might be observable in edge-on spiral galaxies as a 
faint symmetric thickening of the disc due entirely to the disc arms. 

In section (4) we construct an example of the scale-free, isothermal, disc-halo spiral 
system. The initial arms are rigorously discrete in the case of constant gas density, but 
the more likely case has the gas surface density reacting strongly to the spiral potential. 
Examples of this latter variation have been given in figure ([3]). To establish full consistency 
of the model, the magnetohydrodynamics of the gas should be studied in the presence 
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of the system potential. This is rather a comphcated proposition as it is hkely that the 
isothermahty is broken in reahty by the physics of heating and coohng. Moreover the 
sources of the magnetic field are uncertain. 

There is much left undone at this stage even for the isothermal class. This includes 
discussing the possible origin of isothermahty during galactic formation and evolution. This 
is equivalent to asking for the origin of the a = 1 self-similarity, which does seem to arise 
naturally in certain regions of simulated dark matter halos. However this theory and the 
recent simulations do seem to agree on a new picture for at least repeatedly excited spiral 
arms. 
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